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. We discuss a model for the energy distribution and the early space-time evolution of a 
heavy ion collision. We estimate the gluon field generated in the wake of hard processes 
^ ! and through primordial fluctuations of the color charges in the nuclei. Without specifying 
the dynamical mechanism of thermalization we calculate the energy momentum tensor of 
the following plasma phase. The results of this model can be used as initial conditions 
for a further hydrodynamic evolution. 
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1. INTRODUCTION 



Our knowledge about the early phase of a heavy ion collision, is still incomplete. Hydro- 
dynamic calculations compared with data collected at the Relativistic Heavy Ion Collider 
r-| ! (RHIC) suggest a rather early thermalization time Tf < 1 fm/c, at which a deconfined 
quark gluon plasma is created [ . After the time Tf ideal hydrodynamics delivers a 
successful description of the dynamics of the fireball. The initial conditions for the hydro- 
dynamic evolution however remain obscure and are often subject to an educated guess. 

Here we propose a simple model for the initial state that focuses on two important 
questions: how much energy is deposited in the fireball between the receding nuclei, 
and how is this energy distributed in space-time at an early time rj? The model uses 
three simple ingredients: an initial distribution of color sources in the colliding nuclei, 
(ii) the evolution of a classical McLerran-Venugopalan gluon field, (iii) energy-momentum 
conservation. 



2. OUTLINE OF THE MODEL 

Let us choose a momentum scale Qo (which should be much larger than Aqcd)- We 
divide gluon production into hard modes (hard processes with transverse momentum 
Pt > Qo) and soft modes (modes of the gluon field with px < Qo)- Hard modes can 
be rigorously treated within perturbative QCD (pQCD) and implemented in a parton 
cascade. For the soft modes we take guidance from the McLerran-Venugopalan (MV) 
model and describe them by classical solutions to the Yang-Mills equation [D^, F'^'^] = J'^ 
and the continuity equation [D^, J''] = [^. In the MV model the field is generated by 
local color fluctuations in the nuclei which are otherwise color neutral. 

In our case the classical SU (3) current J'^ contains these primordial color distributions 
pi and p2- They are infinitely thin sheets of charge propagating with the speed of light 
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along the light cone from t = — oo to t = 0. Hard processes take place on very short time 
scales T ~ 1/pt- They happen almost instantaneous for the soft fields. For consistency 
color charges induced by hard processes have to be taken into account as additional sources 
of the soft gluon field, which are switched on at t = 0. This is similar to a color ionization 
process: a parton with a certain color is participating in a hard process, leaving a hole 
with the corresponding anticolor. The parton can instead be found in the forward light 
cone. Therefore contains modified color distributions p'l and propagating on the 
light cone from t = to t = +oo and a charge pa in the forward light cone. 

Here we shall not discuss the most general case, but demonstrate the basic features of 
the model for the simple case of a boost invariant current. In an axial gauge, defined 
by the condition A~ + x~ = for the gauge potential A, we consider the current 
Ji(x) = = 1,2) and 

J^{x) = 6{x^) [pi,2(f±) ± e(x±)p(fx)] , (1) 

with p[ 2 = Pi,2 =t P, see Fig. ^ We solve the Yang-Mills equations for small times r. 
Causality ensures that only charges with a transverse distance less than cr can interact 
during this time. This suggests to discretize the charge distributions in the transverse 
plane. We use a rectangular grid of size a for both nuclei. In each bin the classical 
charge with an effective size 1/p < a (p ~ Qo) is put at the center. For ct < a only the 
interaction of charges in the two nuclei which are opposite to each other (i.e. go through 
each other at the time of overlap, t = 0) have to be taken into account. The fireball 
between the receding nuclei is then described by a set of parallel tubes of width a with 
the bins of color charges at both ends (Fig. 12)). 




clei and the tube connecting them which con- 
tains the gluon field generated by the interac- 
Figure 1. The propagation of the color tion. The partons in both bins are represented 
charges defined above on the lightcone. by effective color charges Ti and T2. 



In each bin in the nuclei the color charges of the partons are described by an effective 
classical color charge. It was shown that this classical charge is given by a random walk 
in SU{3) space [0]. The probability distribution for the classical charge T is given by 

w^iim = e-^^\^\'/^ (2) 
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where the length N of the random walk is determined by the number of quarks, antiquarks 
and gluons as = A^^ + Ng + {CA/CF)Ng. 

The Yang-Mills equations can now be solved for each tube using classical charges Ti 
and T2 as sources in the bins on either side. We solve them to all orders in the coupling 
constant g using an expansion in powers of r/i. At small times the dominating field 
component is the longitudinal chromoelectric field F^~ in the tube. It has a finite value 
for r = 0, while the transverse chromoelectric and magnetic fields F*^ vanish initially 
and then increase linearly at early times. The longitudinal chromomagnetic field F*-^ is 
zero at least up to order r^. We have solved the field equations up to order for the 
longitudinal fields and up to order r'^ for the transverse fields. 

Thus at early times the field between the receding nuclei is analogous to a capacitor as 
discussed in The longitudinal chromoelectric field will slow down the charges. The 
trajectory of each bin of partons for t > can be calculated to lowest order by assuming 
that it initially travels with the rapidity Uq of the beam. We can follow the trajectory 
to where it intersects the hypersurface r = rj at space-time rapidities 17* and 772- These 
values limit the longitudinal extent of the tube with the gluon field or plasma (see Fig. 

EI). 

At times cr ~ a the fields of neighboring tubes start to overlap. We add them coher- 
ently, and calculate the resulting energy-momentum tensor T^^ from the gluon field. It 
turns out that the energy-momentum tensor can be parameterized in the following way: 

/ A + C cosh 2?7 Bi cosh rj B2 cosh t] C sinh 2rj \ 

Bi cosh f] A + D E Bi sinh rj 

B2 cosh 7] E A — D B2 sinh t] 

y Csinh2?7 i?isinh?7 B2 sinh rj —A + Ccosh.2r]J 

where A, B = {Bi, B2), C, D and E are functions of time r and the transverse coordinate 
x± but independent of space-time rapidity rj. At early times they behave like A ~ r°, 
B and C ^ respectively, while the functions D and E will not play any further role. 
Using the expansion in r for the gluon field as above we can ensure that d^Tj^ = 0+O{t'^). 

We can now match this field tensor with the tensor of the plasma phase, T^^^ = (e -|- 
p)u^u'^ — pg^^ characterized by a local energy density e, pressure p and a 4-velocity u = 
7(1, v±, Vl) with longitudinal and transverse flow vl and v_i. Without detailed knowledge 
of the decay and thermalization of the gluon field we can still make use of the fact that 
dfjT^^ = where T^,y is the total energy momentum tensor of the fireball which is equal 
to Tj^'^ up to some time before the thermalization time Tf and equal to T^^ after that. 

Here we discuss the case of an instantaneous matching for which T'^^ = Q{tq — t)T^'^ + 
0(r — tq)T!^^ . Energy-momentum conservation implies a set of equations for the plasma: 

"^^Jh^jrlr^' = .+v=(A+c) (i- (^^p^p). ) (4) 

These equations fix transverse flow to be in the direction oi B oc T^*. We also note that 
the longitudinal flow leads to a boost-invariant scenario where rapidity y is equal to 77 up 
to the limiting values ri^2- Of course this is a consequence of the boost-invariant charge 
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distribution in this example. We need another equation to unambiguously determine 
the magnitude of v±, p and e. This can be the equation of state which is used for the 
hydrodynamic evolution for r > Tf. E.g. an ideal gas equation of state would lead to 
e = -{A + C) + V4(A + C)2-3B2. 
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Figure 3. The longitudinal structure of a 
tube in the fireball. The partons in the 
bins to the left and right are decelerated. 
They reach final rapidities ri[ and 773 at 
the thermalization time rj. 
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Figure 4. A sketch of the 2-vector B in 
the transverse plane which determines the 
radial flow in the plasma phase. 



3. SUMMARY 

To summarize, we introduce a model for the early phase of a heavy ion collision that 
can be used to estimate how much energy is deposited by the production of gluon fields 
and how the energy is distributed once the system thermalizes. In this model it is natural 
to have transverse flow and also elliptic flow already at the beginning of the plasma phase. 
Flow originates from the transverse energy flow T^* of the early gluon phase, see Fig. EJ 
Preliminary results have been obtained using a boost invariant charge distribution in the 
nuclei. Some of the details remain to be carried out in more detail. In particular we aim 
for a consistent implementation of hard processes which minimizes the dependence of the 
final results on the hard scale Qo- Hydrodynamic simulations using our results are in 
preparation. 
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